Phase transition for the Ising model on the 
Critical Lorentzian triangulation 

Maxim Krikun*^'^ Anatoly Yambartsev''"'^ 
October 13, 2008 



1 Institut Elie Cartan Nancy (lECN), Nancy-Universite, CNRS, INRIA, Boulevard des 

Aiguillcttes B.P. 239 F-54506 Vandoeuvre-les-Nancy, France. 

E-mail: krikun@iecn.u-nancy.fr 

current aiEliation: Google Ireland, Dublin, Ireland 

^ Department of Statistics, Institute of Mathematics and Statistics, University of Sao 
Paulo, Rua do Matao 1010, CEP 05508-090, Sao Paulo SP, Brazil. 
E-mail: yambar@ime.usp.br 

Abstract 

Ising model without external field on an infinite Lorentzian tri- 
angulation sampled from the uniform distribution is considered. We 
prove uniqueness of the Gibbs measure in the high temperature region 
and coexistence of at least two Gibbs measures at low temperature. 
The proofs are based on the disagreement percolation method and on 
a variant of Peierls method. The critical temperature is shown to be 
constant a.s. 

Keywords: Lorentzian triangulation, Ising model, dynamical triangu- 
lation, quantum gravity 

AMS 2000 Subject Classifications: 82B20, 82B26, 60J80 



* Current affiliation: Google Ireland, Dublin, Ireland. 
tPartially supported by BQR 2007 UMR7502 lAEM 0039. 

tPartially supported by the "Rede Matematica Brasil-Franga, CNPq (306092/2007-7) 
and CNPR "Edital Universal 2006" (471925/2006-3). 



1 



1 Introduction 



Triangulations, and planar graphs in general, appear in physics in the context 
of 2-dimensional quantum gravity as a model for the discretized time-space. 
Perhaps the best understood it the model of Euclidean Dynamical Triangu- 
lations, which can be viewed as a way of constructing a random graph by 
gluing together a large number of equilateral triangles in all possible ways, 
with only topological conditions imposed on such gluing. Putting a spin sys- 
tem on such a random graph can be interpreted as a coupling of gravity with 
matter, and was an object of persistent interest in physics since the successful 
application of matrix integral methods to the Ising model on random lattice 
by Kazakov [S]. 

More recently, a model of Casual Dynamical Triangulations was intro- 
duced (see [9] for an overview). The distinguishing feature of this model is 
its lack of isotropy — the triangulation now has a distinguished time-like 
direction, giving it a partial order structure similar to Minkowski space, and 
imposing some non-topological restrictions on the way elementary triangles 
are glued. This last fact destroys the connection between the model and ma- 
trix integrals, in particular the analysis of the Ising model requires completely 
different methods (see e.g. [3]). 

From a mathematical perspective, we deal here with nothing but a spin 
system on a random graph. Random graphs, arising from the CDT approach, 
were considered in [llj under the name of Lorentzian models. In the present 
paper we consider the Ising model on such graphs. When defining the model 
we pursue the formal Gibbsian approach [5]; namely, given a realization of 
an infinite triangulation, we consider probability measures on the set of spin 
configurations that correspond to a certain formal Hamiltonian. 

Our setting is drastically different from e.g. [8] and [3J in that we do not 
consider "simultaneous randomness" , when both the triangulations and spin 
configurations are included into one Hamiltonian. Instead we first sample 
an infinite triangulation from some natural "uniform" measure, and then 
run an Ising model on it, thus the resulting semi-direct product measure is 
"quenched" . 

A modest goal of this paper is to estabhsh a phase transition for the Ising 
model in the above described "quenched" setting (the "annealed" version of 
the problem is surely interesting, but is also more technically challenging, so 
we don't attempt it for the moment). In Section 3 we use a variant of Peierls 
method to prove non-uniqueness of the Gibbs measure at low temperature. 
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Quite surprisingly, proving the uniqueness at high temperature is not easy 
- the difficulty consists in presence of vertices of arbitrarily large degree, 
which does not allow for immediate application of uniqueness criteria such 
as e.g. [H]. We resort instead to the method of disagreement percolation |T3], 
and use the idea of "ungluing", borrowed from the paper [2], to get rid of 
vertices of very high degree. 

Finally, in Section 5 we show that the critical temperature is in fact 
non-random and coincides for a.e. random Lorentzian triangulation. Sec- 
tion 2 below contains the main definitions and summarizes some of the results 

aim- 

We thank E. Pechersky for numerous useful discussions during the prepa- 
ration of this paper. 

2 Definitions and Main Results 

Now we define rooted infinite Lorentzian triangulations in a cylinder C = 
X [0,oo). 

Definition 2.1. Consider a connected graph G embedded in a cylinder C. A 
face is a connected component ofC\G. The face is a triangle if its boundary 
meets precisely three edges of the graph. An embedded triangulation T is 
such a graph G together with a subset of the triangular faces of G. Let 
the support S{T) d C be the union of G and the triangular faces in T. 
Two embedded triangulations T and T' are considered equivalent if there is a 
homeomorphism of S{T) and S{T') that corresponds T and T' . 

For convenience, we usually abbreviate "equivalence class of embedded 
triangulations" to "triangulation". This should not cause much confusion. 
We suppose that the number of the vertices of G is finite or countable. 

Definition 2.2. A triangulation T of C is called Lorentzian if the following 
conditions hold: each triangular face of T belongs to some strip x [j', j + 
= 0,1,..., and has all vertices and exactly one edge on the boundary 
{S^ X {]}) U {S^^ X {i + 1}) of the strip x + 1]; and the number of 
edges on x {j} is positive and finite for any j = 0, 1, . . . . 

In this paper we will consider only the case when the number of edges 
on the first level x {0} equal to 1. This is not restriction, only it gives 
formulas more clean. 
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Definition 2.3. A triangulation T is called rooted if it has a root. The root 
in the triangulation T consists of a triangle t ofT, called the root face, with 
an ordering on its vertices (x, z) . The vertex x is the root vertex and the 
directed edge {x,y) is the root edge. The x and {x,y) belong to x {0}. 

Note that this definition also means that the homeomorphism in the def- 
inition of the equivalence class respects the root vertex and the root edge. 
For convenience, we usually abbreviate "equivalence class of embedded rooted 
Lorentzian triangulations" to "Lorentzian triangulation" or LT. 

In the same way we also can define a Lorentzian triangulation of a cylinder 
Cat = 5*^ X [0, N]. Let LT^r and LTqo denote the set of Lorentzian triangula- 
tions with support Cat and C correspondingly. 

2.1 Gibbs and Uniform Lorentzian triangulations. 

Let LTat be the set of all Lorentzian triangulations with only one (rooted) 
edge on the root boundary and with slices. The number of edges on the 
upper boundary 5*^ x {A^} is not fixed. Introduce a Gibbs measure on the 
(countable) set hTjy: 

P^,^(T) = Z[o,V]exp(-/iF(T)), (2.1) 

where F(T) denotes the number of triangles in a triangulation T and Z^Q^j^j 
is the partition function: 

The measure on the set of infinite triangulations LTqo is then defined as a 
weak limit 

:= lim Pat,^ . 

It was shown in that this limit exists for all fi > ficr '■= In 2. 

Theorem 1 (IH]). Let kn be the number of vertices at n-th level in a trian- 
gulation T for each n > 0. 

• For /i > ln2 under the limiting measure P^ the sequence {kn} is a 
positive recurrent Markov chain. 



4 



• For fi = = ln2 the sequence {A;„} is distributed as the branching 
process ^„ with geometric offspring distribution with parameter 1/2, 
conditioned to non- extinction at infinity. 

Below we briefly sketch the proof of the second part of Theorem [H a 
deeper investigation of related ideas will appear in [12]. 

Given a triangulation T G LT^r, define the subgraph r C T by taking, 
for each vertex u G T, the leftmost edge going from v downwards (see fig. [I]). 
The graph thus obtained is a spanning forest of T, and moreover, if one 
associates with each vertex of r it's height in T then T can be completely 
reconstructed knowing r. We call r the tree parametrization of T. 




Figure 1: Tree parametrization 



For every vertex u G t denote by 6^ it's out-degree, i.e. the number of 
edges of r going from u upwards. Comparing the out-degrees in r to the 
number of vertical edges in T, and comparing the latter to the total number 
of triangles, it is not hard to obtain the identity 

J2 {Su + 1) = F{T) (2.2) 

«eT\5ix{Af} 

(the sum on the left runs over all vertices of r except for the A^-th level). 
Thus under the measure PN,iJ.cr the probability of a forest r is proportional 

to 

n (^)''" (2-3) 

ueT\sT-x{N} 

which is exactly the probability to observe r as a realization of a branching 
process with offspring distribution Geom(l/2). After normalization we'll 
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obtain, on the left in fl2.3l) . the probabihty PN,^i{'T) as defined by fl2.1l) . an on 
the right the conditional probabihty to see r as a reahzation of the branching 
process ^ given > 0. So quite naturaUy when N —>■ oo the distribution 
of r converges to the Gahon- Watson tree, conditioned to non-extinction at 
infinity. 

In particular it follows from Theorem [T] that 

P^AK = m) = Pr{^^ = m I > 0) = mPr(a = m). (2.4) 

Remark 2.1. The last equality in fl2.4p means that the measure P^^,, on tri- 
angulations can be considered as a Q-process defined by Athreya and Ney 
yy for a critical Galton- Watson branching process. Such a process is exactly 
a critical Galton- Watson tree conditioned to survive forever. 

We will also note for further use that the offspring generating function of 
the branching process ^, 

fc>0 

and the generating function for C,n (with initial condition = 1) is 

n-{n-l)s 

w) = — —. • 

n + 1 — ns 

2.2 Ising model on Uniform Infinite Lorentzian trian- 
gulation — quenched case. 

Let T be some fixed Lorentzian triangulation, T G LTqo. Let T/v be the 
projection of T on the cylinder Cat. We associate with every vertex v a spin 
(Jv £ !}• Let S(T) and Sjv(T) denote the set of of all spin configurations 
on T and Tjv, respectively. The Ising model on T is defined by a formal 
Hamiltonian 

H{a)= J2 (2-5) 

{v,v')ev 

where {v, v') means that vertices w, v' are neighbors, i.e. are connected by 
an edge in T. Let OTn be the set of vertices of T that lie on the circle 
5*^ X {A^ -|- 1}. Fix some configuration on the boundary STat and denote it 
da. The Gibbs distribution with boundary condition da is defined by the 
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following. Let V{Ti^) be the set of all vertices in T^, then the energy of 
configuration a G T,n(T) is 



(2,6) 



{v,v'): v,v'&ViTN) 



(v,v'}:veViTN),v'£aTN 



which defines the probability 




exp { — /^i^Ar (a I (9(j) } 

ZN,da{T) 



(2.7) 



where 




J2 exp{-pHr,{a\da)}. 



<TeSjv{T) 

When N ^ oo, for any sequence of boundary conditions da, a limit (at least 
along some subsequence) of measures Pj^ g„ exists by compactness. Such a 
limit is a probability measure on S(T) with a natural a- algebra, which we 
refer to Gibbs measure. 

In general, it is well known that at least one Gibbs measure exists for the 
Ising model on any locally finite graph and for any value of the parameter 
P (see, e.g., [6] page 71). It is also known that the existence of more than 
one Gibbs measure is increasing in j3, i.e. there exists a critical value jSc G 
[0, oo] such that there is a unique Gibbs measure when (3 > jSc, and multiple 
Gibbs measures when (3 < (3c (see [7j for an overview of relations between 
percolation and Ising model on general graphs). 

Thus when considering the Ising model on Lorentzian triangulations it is 
natural to ask whether the critical temperature is finite (different from both 
and oo), and whether it depends on the triangulation. In the following two 
sections we show that the critical temperature is a.s. bounded both from 
and oo. In the last section we prove that the critical temperature obeys a 
zero-one law and is therefore a.s. constant. 

3 Phase transition at low temperature 

In this section we prove the following theorem. 

Theorem 2. There exists a (3q such that for all (3 G {(3q,oo) there exist at 
least two Gibbs measures for P,, -a.e. T. 
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We remind first the classical Peierls method for the Ising model on . 
Let A C be a large square box centered at the origin, and let be the 
distribution for the spins in A under the condition that all of the spins outside 
of A are negative. When the size of the box tends to infinity, converges 
to some probability measure P~, which we call Gibbs measure with negative 
boundary conditions. The measure is defined similarly taking positive 
boundary conditions. 

Due to an obvious symmetry between P^ and P^ we have 

p-(ao = +l)=P+((To = -l), 

therefore if the Gibbs measure is unique then necessarily P^ = P^ and 

p-((To = +l)=P'(ao = -l) = l/2. (3.1) 

Thus our goal will be to disprove (13.11) for large enough j3. 

Consider some configuration a such that ctq = +1 and cTj = — 1 for all 
i ^ A. Define the cluster C as the maximal connected subgraph of Z^, 
containing the origin, such that cXj = +1 for all i G K+. Let 7 be the contour 
in the dual graph (Z^)*, corresponding to the outer boundary of K^, and let 
a' be the configuration obtained from a by inverting all the spins inside 7. 
Note that ctq = —1. Let C be the set of all contours in (Z^)* surrounding the 
origin, and consider the apphcation 

inv : S|o-o=+i S|o-o=-i x C , s 

Clearly inv is injective: given {a', 7) one can easily reconstruct a by inverting 
the spins inside 7. Also the following property holds: 

Pj,{cr) = P^{a')-e-'^\^\, (3.3) 

where I7I denotes the length of 7. Indeed, for every edge traversed 
by 7 we have ai 7^ aj by construction, but after inversion a- = cr'j, so the 
contribution of the pair {i,j) to the Hamiltonian is increased by 2. On the 
other hand, for all other pairs the contribution to the Hamiltonian 

doesn't change. Informally we can write H{a) — H{a') = 2\'~f\, which is 
equivalent to (13. 3p . 
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Now taking a sum over all cr : o-q = +1 in fl3.3l) we get 



(t:(to=+1 



(3.4) 



where C„ is the set of contours 7 G C of length exactly n. 

But < n3": indeed, every contour 7 G C„ intersects the x-axis 

somewhere between and n. Starting from this intersection point there are 
at most 3" distinct self-avoiding paths of length n in (Z^)*, only a few of 
them really belonging to C„. Therefore we have the inequality 



PA(ao = +l)<PA(ao = -l)-5^e- 



(3.5) 



n>l 



and by taking f3 large enough the sum in the right-hand side of (13. 5p can be 
made strictly less than one. Taking the limit A ^ we establish 



-1) < P"(ao 



-1). 



Thus P 7^ P+ and the Gibbs measure for the Ising model in at inverse 
temperature j3 is not unique. 





Figure 2: Peierls' method on Z^ (left) and on a random Lorentzian triangu- 
lation (right). Dashed line represents the inversion contour 7. 



Next we will slightly generalize the above classical argument. 
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Lemma 3.1. Let G be an infinite planar graph and G* it's planar dual. Let 
Vq he a vertex of G and denote by C„(fo) be the set of contours of length n 
in G* , separating vq from the infinite part of the graph. 
If the following sum is finite, 

$^#C„(t;o)e-'^"<oo, (3.6) 

n>l 

then the Gibbs measure for the Ising model on G at inverse temperature j3 is 
not unique. 

Remark 3.1. Strictly speaking, in the above Lemma we should also require 
the graph G to be one-ended; i.e. for any finite subgraph B G G the com- 
plement G\B must contain exactly one infinite connected component. The 
reason for this is that if the graph G fails to be one-ended (e.g. a doubly- 
infinite cylinder) then the outer boundary of the cluster may happen to 
have multiple connected components, which slightly complicates the argu- 
ment. We don't consider this case in detail since Lorentzian triangulations 
are one-ended by construction. 

Proof. Assume that (13. 6p holds. Then for some large 

#C„(t;o)e-'^" < 1. 

n>N 

Also, for every n the set C„(fo) is finite, so there exists some large R such 
that the ball Bji{vq) contains all of the Gn{vo), n = 1 . . . N. 

Let now and be the Gibbs measures for the Ising model on G, 
constructed with positive and negative boundary conditions respectively, and 
let us compare the events 

A+ := (cr^ = +1 for all v e Br{vo)) 

and 

A_ := {ay = -1 for all v G Br{vo)). 
Proceeding as in the case of above, we can show that 

P-(A+) < p-iA.) • 5^e-2^"#C„(i?H(t;o)), 

n>l 
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where C„(i?/j(fo)) now is the set of contours of length n in G* that surround 
the whole ball Bji{vq). But by construction we have Cn{Bfj{vQ)) = for 
n < N and C„(-Bk(^o)) C. C„(fo) for all other n, therefore 

^e-2/^"#C„(5^(t;o)) < 5^ e-2/5"#C„(t;o) < 1 

n>l n>N 

and 

P-(A+)<P-(A_). 
Since by symmetry = P"(A+), it follows that P" ^ F+. □ 

Lemma 3.2. Let T be a random Lorentzian triangulation, and let Vq he the 
root vertex ofT. There exists such that for every (3 > [3q 

$^#C„(t;o)e-'^"<oo P^,,-a.e, (3.7) 

n>l 

where i^Cn{vo) is defined as in Lemma \3. 1[ 

Proof. In order to prove (13.71) it will be sufficient to show that the expectation 
(with respect to the measure P^^^) of the sum is finite 

oo 
n=0 

First of all we choose in any Lorentzian triangulation T a vertical path 
7oo = 7oo(^) = (^^0) Vi, . . .) starting at the root Vq and such that Vi E S^x {i}. 
Let Cji^ni^o) C C„(fo) be the set of contours of length n which surround Vq 
and intersect 7oo at height R. Note that any such contour does not exit from 
the strip [R — n, R + n]. Let also Sr^u be the number of particles in the tree 
parametrization of T at height R — n which have nonempty offspring in the 
generation located at height R + n. Since every contour from Cr^ji, in order 
to surround vq, must cross each of the Sr^u corresponding subtrees, we have 

{SR^n >n) ^ {CR^n = 0). 

On the other hand 

since the contours C/j^„ live on the dual graph T*, which has all vertices of 
degree 3, thus there are at most 2^ self- avoiding paths with a fixed starting 
point (which is in our case the intersection with 7oo)- 
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Therefore we have 

oo oo oo 

E... *CMe-'^- < e-'''"2" ^ P.^Sn^n < n) (3.9) 

n=0 n=0 R=l 

and for this sum to be finite it's enough to show that sum over R has poly- 
nomial order in n. First let us estimate 

oo oo 

Y ^^^r{SR,n < n) < + ^ P,.XSR,n < n). 

R=l R>rr> 

and write for some e G (0, 1/4) 

P^l,XSR,n <n) = Y ^t^oX^R^n < n I kR_n = k)Pf,^^{kR_n = k) 

+ Y ^^^cr{SR,n < n I kR.n = k)P^^^{kR.n = k%3.10) 

k>R'^ 

The first sum above satisfies the inequality 

Y ^^^rikR-n = k) = P,AkR.-n < R') < {j^f ■ (S-H) 

k<R<^ 

Indeed, thanks the representation (12.41) the generating function ipn{s) of kn is 
related to the generating function ipn{s) of the branching process ^„ described 
above by the following equation 

¥.„(s) = #;(s). (3.12) 

Since ipn{s) can be calculated explicitly 

n — (n — l)s 

Ms) = 3.13 

n + I — ns 



from (I3.12P and (I3.13P we obtain 



^"(^) = 7TT ^' ^3.14) 



P,Jk^ = k) = ^^^n(.)^^^ 

kn^^^ k 



(n + l)'^ ' 
12 



< ^, (3.15) 



k ^ 

which proves the estimate (13. lip . Since e < 1/4 there exists a constant ci 
such that 

The second sum in (13.101) can be estimated by the probabihty P^^^{SR^n < n \ 
kR-n = -R*^)- To do this, let us remind some properties of conditioned or size- 
biased Galton- Watson trees (see e.g. [10], [1] or [1]). If r is a Galton- Watson 
tree of a critical branching process, conditioned to non-extinction at infinity, 
then conditionally on kR^n the subtrees of r, originating on the {R — n)-th 
level, are distributed as follows: one subtree, chosen uniformly at random, has 
the same distribution as r (i.e. is infinite), while the remaining subtrees are 
regular Galton- Watson trees corresponding to the original branching process 
(and are finite a.s.). 

The probability for a particle of the branching process ^ to survive up to 
time 2n equals l/{2n+ 1), thus conditionally on /cR-n, Sr^u is stochastically 
minorated by the binomial distribution with parameters {kR-n, 2^;^)- Using 
Hoeffding's inequality for binomial distribution we obtain 

P,ASR,n < n I kR.^ = < exp(-2 ^^^/^^^+/^~^^' ) (3.17) 

and the sum over R> is bounded by some absolute constant C2 

/ (i?7(2n + 1) -n)2\ 2^ / 2R^ \ , , 

Thus we have proved that 

CO oo 

^ 5^e-2/3"2"(n5 + ci + c^) (3.19) 

n=0 n=0 

and so for all /? > ^ the sum (13. 7p is finite. □ 
Proof of Theorem [H Theorem [2] follows immediately from Lemmas 13.11 
and [321 
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4 Uniqueness for high temperature 



In this section we prove the following theorem 

Theorem 3. There exists a small enough (3^ such that for every (3 G [O,/?/^) 
for P^^^-a.e. Lorentzian triangulation T the Gibbs measure for the Ising model 
on T is unique. 

The proof goes as follows. First, for a fixed triangulation T we use dis- 
agreement percolation to reduce the problem of uniqueness of a Gibbs mea- 
sure to the problem of existence of an infinite open cluster under some specific 
site-percolation model on T. Then we consider the joint distribution of a tri- 
angulation with this percolation model on it (i.e. with randomly open/closed 
vertices), where the triangulation T is distributed according to P^^^, and con- 
ditionally on T the state of each vertex in T is chosen independently (but 
with vertex-dependent probabilities). We show that the probability for such 
a random triangulation to contain an infinite open cluster is 0, therefore this 
probability is for P^^^-a.e. triangulation as well, and the theorem follows. 

4.1 Disagreement percolation 

The following result of Van den Berg and Maes [13] provides a sufficient 
condition for the uniqueness of the Gibbs measure, expressed in terms of 
some percolation-type problem. 

Theorem 4. Let G be a countably infinite locally finite graph with the vertex 
set V, and let S be a finite spin space. Let Y be a specification of a Markov 
field on G, i.e. for every finite A C V, YA{-,ri) is the conditional distribu- 
tion of the spins on A, given that the spins outside of A coincide with the 
configuration rj G . 

Define for every v E V 

py := max pvar(yM(-, r]), Y|^}(-, V))> 

where p^^^y- denotes the distance in variation. Let P^'^ be the product mea- 
sure under which each vertex v is open with probability p^ and closed with 
probability 1 — p^. 

If under site-percolation on G, one has 

P^'^ {there is an infinite open path} = 0, 

then the specification Y admits at most one Gibbs measure. 
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In the case of the Ising model, the only spins that affect the distribution 
of the spin at a given vertex v are those at vertices u connected to v by an 
edge (denoted u ^ v). Let be the degree of v and let = J2ur^v^u- The 
conditional distribution of ay, given the values of spins cr„ for all m ~ f , only 
depends on Sy and is given by 



P^{ay = +l\Sy} 



P'^iay = -l\Sy} 



The distance in variation between two such conditional distributions is max- 
imized by taking two extreme values of Sy, namely Sy = dy and Sy = —dy, 
corresponding to all-plus and all-minus boundary conditions. Therefore we 
have 



Pv 



1 

2 

+ 



'iv _j_ 

o-Pdv 



g Pd^ _j_ g/: 
^(3dv 



Q—Pdv _|_ g/Sdi, 

and we need to prove the following 



g-/3rf. + g/S 



tanh(/3(i^) 



Lemma 4.1. Let T he an infinite random Lorentzian triangulation sampled 
from the measure P and let each vertex v E T be open with probability 
Py = ta.iah{Pdy) and closed otherwise. Then for small enough [3 the probability 
that there is an infinite open path in T is zero. 



4.2 Elementary approach to non-percolation 

Imagine for a moment that the graph G that we deal with in Lemma 14.11 
has uniformly bounded degrees. If with positive probability there exists an 
infinite open percolation cluster, then with positive probability this cluster 
will also be connected to the origin. Consider the event 

^ r there exists a self-avoiding open path of length n 1 
" \ starting at the origin. J 

The number of paths of length n, starting at the origin, is at most D", 
where D is the maximal degree of a vertex in our graph. Obviously the same 
estimate is valid for the number of self-avoiding paths. The probability for a 
vertex to be open is bounded by tanh(/3D) < (3D, therefore we can estimate 

P'^'^i.An) < L'"tanh(/5L')" < /^"L'^". 
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Taking (3 < we get PP'^(A„) ^ as 72 oo, and therefore the 

probabihty that the origin belongs to an infinite open cluster is zero. 

Unfortunately, the above argument can't be applied directly in the con- 
text of Lemma HTTl in an infinite random Lorentzian triangulation with prob- 
ability one one will encounter regions consisting of vertices of arbitrarily large 
degrees, and these regions can be arbitrarily large as well. On the other hand, 
we know that in a Lorentzian triangulation the degree of a typical vertex is 
not too large (in some sense, which we're not going to make precise in this 
paper, it can be approximated by a sum of two independent geometric dis- 
tributions plus a constant), therefore one can hope to show that the vertices 
of very large degree are so rare that they don't mess up the picture. 

Before proceeding further, let us remark that in the above argument one 
can consider, instead of the set of all the self-avoiding paths, only the paths 
that are both self-avoiding and locally geodesic. We say that a path 7 in 
a graph G is locally geodesic if whenever two vertices, belonging to 7, are 
connected by and edge in G, this edge also belongs to 7 (in other words, a 
locally geodesic path avoids making unnecessary detours). Clearly, if under 
site percolation on G there exists an infinite open self-avoiding path starting 
at vq, there exists also an infinite self- avoiding locally geodesic path. It is 
also easy to see that if 7 is a locally geodesic path, then every vertex t> G 7 
has at most two neighbors belonging to 7. We will use this observation in 
the next section. 

4.3 Generalization to random triangulations 

Fix an infinite Lorentzian triangulation T and let 7 be a finite self-avoiding, 
locally geodesic path in T, starting at the root vertex (which we denote by 
Vo). Let r C T be a 1-neighhorhood oi 7, i.e. a sub-triangulation consisting of 
the path 7 together with all triangles, adjacent to 7. Such a sub-triangulation 
enjoys the rigidity property: if F can be embedded into some Lorentzian tri- 
angulation, then this embedding is uniqu^. We define the event "T contains 
F" (denoted F C T) to be the subset of triangulations in LTqo for which such 
an embedding exists. 

Lemma 4.2. Let T he a 1 -neighborhood of a locally geodesic path 7, let 7 
have length n and let the vertices of 'y have degrees d^, di, . . . , d^. 

^ we require that the embedding maps vq to the root vertex of T, and sends horizon- 
tal/vertical edges of T to edges of the same type in T 
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The probability that a random Lorentzian triangulation T contains T is 
bounded by 

n 
i=0 

for some absolute constant C. 

We can use this lemma to obtain estimates on percolation probability as 
follows. Given the degrees d^, . . . ,dn one can completely specify F using the 
following collection of numbers 

• for each vertex we need to specify how it's degree dj is split into d^"^^ 

and d^^"'\ i.e. edges going up and down from the vertex — this gives 
one integer parameter in the range [1, dj]; 

• also we need to specify the edge along which the path 7 entered the 
vertex and the edge along which it left; this adds two more parameters. 

The total number of meaningful combinations for the above parameters 
clearly does not exceed Yl^^- Now, if T is a P^^^-distributed random tri- 
angulation with site-percolation on it, by Lemma [4.21 the probability for T 
to contain an open path from the origin with vertex degrees do, . . . ,dn is 
bounded by 




Summing over all possible values oi do, . . . ,dn we obtain the estimate 

do,...,dn ■> d>i 

for the probability to have an open path of length n in a random triangula- 
tion. For (3 sufficiently small, the last expression tends to as n — > 00, and 
the probability to have an infinite open path from the origin is zero. 

Note that here we are estimating the annealed probability (semi-direct 
product P^^^ xPP''^) of the existence of open path. Clearly, if this annealed 
probability is zero, then also for P^^^-a.e. triangulation T the P^'^-probability 
of an infinite open path is zero as well. 

Therefore it remains to prove Lemma [4. 2[ 
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4.4 Elementary perturbations 



Consider a finite Lorentzian triangulation T G LTat containing an internal 
vertex v of very liigfi degree, say > 100. Let dl"^^ and di^^^ be tlie degree 
contributed to dy by edges going upwards and downwards from v, so that 
4 = 4"P) + rfl''")+2. LetO{v) be a 1-neighborhood of f , i.e. all the triangles 
adjacent to v, and let T' be the complement to 0{v) in T. 

Let us estimate P n,ij.{T\T') , i.e. the probability to see T as a realization 
of a random Lorentzian triangulation, conditionally on the event that this 
triangulation contains T'. For this purpose consider a modification of the 
1-neighborhood of v, which consists in adding a pair of triangles as shown on 
fig. [31 There are d^v^^^ x di^"^ to distinct ways to make such a modification. 




Figure 3: Inserting a pair of triangles at a high-degree vertex v. The position 
to insert new triangles is determined by the choice of two neighboring vertices 
on the layers below and above v. 



and every time the Hamiltonian of the Gibbs measure (12. ip is increased by 
2/i (since there are 2 triangles added) 
Therefore we can estimate 



PnAT I T') < -T—r-TTT < 



2m g2M 



A better estimate can be obtained by applying k elementary perturbations 
at V simultaneously. Such a set of perturbations is specified by a sequence 
of pairs Mi'^"'')}f=i of vertices lying on the layers above and below v 

(fig-SD- 

In order to assure that all the perturbations can be applied unambigu- 
ously, we require that the vertices u^^'^\ . . -u^^^^ to be ordered from left to 
right, as well as u'i'^\ . . . u^^^\ We may however allow some vertices to coin- 
cide - otherwise it would be impossible to apply n perturbations at once in 
a vertex whose up-degree or down-degree are not large enough. 
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Figure 4: Multiple insertions at the same vertex. Note that one of the edges 
is used twice for insertions, and the corresponding vertex has two dots. 



If is large, one of d^^^ or d^"'^ is at least {d^ — 2)/2, and therefore the 
number of distinct ways to choose, say, k — 10 modifications is at least 



{d, - 2)/2\ 
10 J 



Ci^d^^ + 0{dl). 



Now consider a path 7 C F with vertices of degrees do, • • • , c?n- Let Wy 
be the set of possible modifications we can make to F, by inserting exactly 
10 pairs of triangles at vertices of high degree (say, di > 100), and leaving all 
other vertices intact. Then for some constant Cio 

#Wr>l[Ciodf. 

Let LTr,Ar C LTjv be the set of Lorentzian triangulations that contain F. 
Clearly, each modification w G Wr can be applied to every T G LTr,Ar, 
producing a modified triangulation which we denote by w{T). Consider the 
application of modifications w to a triangulation T as a mapping 

app: LTr,ivxVFr ^ LT^ 

{T,w) ^ w{T). ^^-^^ 

Since any modification to T G LTp^jv adds at most 20n new triangles, we 
have, for all (T, w) G LTr x Wr 

Taking a sum over all (T, w) we get 

J2 PNA^iT)) > e-^'-^^Wr (4-2) 

{T,w)&JTr,N^Wr TeLTr.iv 

= e-2°"'^#iyrPiv,^(r). 
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4.5 Overcounting via random reconstruction 

In order to get an estimate on PAr^^(r) from (14.21) . it would be helpful if the 
application (14. ip was an injection. But this is not easy to prove (if at all 
true!), so instead we will estimate the overcounting 



Lemma 4.3. Given a path 7 with 1-neighborhood T, with vertex degrees 
do, ... , dn, we have for any T' G LTjy 



Proof. First note that each elementary perturbation (i.e. insertion of two 
adjacent triangles) can be undone by collapsing the newly inserted horizontal 
edge. Therefore in order to reconstruct the pair (T, w), knowing T' = w{T), 
it will be sufficient to identify the edges of T' that were added by w. 

It's not clear whether we can perform such a reconstruction determinis- 
tically, but we certainly can achieve it with some probability with the help of 
the following randomized algorithm: 

• the algorithm starts from the root vertex and at each step moves to an 
adjacent vertex, chosen uniformly at random; 

• it also modifies the triangulation along his way as following 

• each time arriving to a new vertex it should decide whether to contract 
a sequence of 10 horizontal edges (at least one of these edges must be 
adjacent to the current location), or do nothing. This gives 12 = 11 + 1 
possibilities, among which one is chosen uniformly at random. 

• after n steps, the path traversed by the algorithm is declared to be the 
(conjectured) path 7, and if it really is, and if the resulting triangu- 
lation coincides with the original triangulation t, the reconstruction is 
considered successful. 

Most of the time this algorithm will fail, but with some probability it may 
actually reconstruct both T and w. Let us now estimate this probability. 

First note that the degree of vq in T' can be larger than the degree of 
Vq in T because of the modifications made to the second vertex of 7, but in 
any case it will not exceed c/q + 10, since each elementary modification in vi 



overr(r') := #{(T, w) G LTr,^ x Wr \ w{T) = T'}. 




(4.3) 
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increases the degree of Vq at most by 1. Therefore, the algorithm will guess 
the correct direction from vq with probability at least l/((io + 10). 

Assuming the direction it took on the first step was correct, there is 
exactly 1 right choice out of 12 on the next step (collapse some 10 nearby 
horizontal edges, which can be chosen in 11 different ways, or do nothing), 
so the probability not to fail is 1/12. 

Assuming again that the modifications to the second vertex vi were un- 
done correctly, there is at least one edge that will lead to the next vertex V2- 
Since the path 7 is locally geodesic, at most two neighbors of vi belong to 
7, namely these are vq and V2- Modifications at V2 add at most 10 to the 
degree of Vi, so with probability at least l/{di + 10) the algorithm will take 
the correct direction again. 

Proceeding this way, we see that the probability to win is at least 

n 

Now for every T' and for every {T,w) G LTr,Ar x Wr such that w{T) = 
T' the probability to correctly reconstruct {T,w) from T' is at least (14. 4p . 
therefore 

overr(r') < JJ 12(rfj + 10). 

This finishes the proof of Lemma 14.31 □ 
Finally, combining the inequalities (14. 2 p and (14. 3 p we get 

e-20"/^#l^r P^,^(r) < J2 ^NA^iT)) 

(T,«))GLTr,jvxH/r 

< Yl overr(T') P;v,M(r') 

n 

< l[l2{d, + 10) PnAT') 

j=0 T'eLT 



< 



and 

n 



J]l2(rf, + 10), 

j=0 



e^o- n 12(rf. + 10) „ 

p-.(r)< — — <n 



j=0 
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Since the last estimate is uniform in N, it holds also for the limiting mea- 
sure P^, and setting fi = fi^r we obtain the statement of Lemma [4.2[ This 
also finishes the proof of Theorem [31 

5 Critical temperature is constant a.s. 

Consider the critical temperature l3c{G) of the Ising model on a graph G as 
a function of G. In the above two sections we have shown that when T is a 
P^^^-random Lorentzian triangulation, we have [3c{T) G [/3/i,/5o] a-S- In this 
section we show that in fact 

Theorem 5. (3c{T) is constant P^^^-a.s. 

The proof relies on two lemmas. 

Lemma 5.1. Let G, G' be two locally finite infinite graphs that differ only 
by a finite subgraph, i.e. there exist finite subgraphs H G G, H' G G' such 
that G\H IS isomorphic to G'\H' . Then f3c{G) = f3c{G'). 

Proof. This statement is widely known and follows from Theorem 7.33 in [5]. 

More exactly, consider the Ising model on G at some fixed inverse tem- 
perature [3 e (0, cxo). Let Y be the corresponding specification (collection 
of conditional distributions) and denote by Q{Y) the set of Gibbs measures, 
satisfying this specification. It is an easy fact that GiX) is a convex set in 
the linear space of finite measures on the set of configurations. 

Any specification F on G can be naturally restricted to a specification Y 
on G\H. Namely, for each finite A G G\H and each r] G {—1, l}*"^^^ put 

Ya{;v) = P{a, = {-)iVieA;ai = r]iVie{G\H)\A}. 

J2 W((-)U^h) (5.1) 

(the last sum runs over all spin configurations in H). Informally, the re- 
stricted specification Y can be interpreted as an Ising model on G, but with 
the spins in H being hidden from the observer. 

Let now Y' be the specification for the Ising model on G', and let Y' be an 
analogous restriction of Y' to G'\H'. Since G\H and G'\H' are isomorphic, 
the specifications Y and Y' are defined over the same graph, and it follows 
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from fl5.ll) that they are equivalent in that there exists a constant c > 1 such 
that 



for all possible values of A, t] and (•). Theorem 7.33 in [6] then implies that 
Q{Y) is affinely isomorphic to Q{Y'). In particular G{Y) consists of a single 
element - a unique Gibbs measure - if and only if Q{Y') does. 

Since any Gibbs measure on spin configurations in G\H, satisfying Y, 
can be uniquely extended to a Gibbs measure on G, satisfying Y, and the 
same holds for C, H', Y', the lemma follows. 

Lemma 5.2. A Galton-Watson tree for a critical branching process, condi- 
tioned to non- extinction at infinity, can he described as following 

• it contains a single infinite path (a spinej, {fo,t'i, . . .}, starting at the 
root vertex; 

• at each vertex Vi of the spine a pair of finite trees {Li,Ri) is attached, 
one of each side of the spine; 

• and the pairs {Li, Ri) are independent and identically distributed. 

Remark 5.1. In fact, the law of {Li,Ri) can be explicitly expressed in terms 
of the original branching process, but for our purposes it's enough to know 
that the trees are i.i.d. 

Proof. See e.g. fTO] or [i]. 

Proof of Theorem [51 Consider now the tree parametrization r of a P^^^- 
distributed infinite Lorentzian triangulation T. By Lemma 15.21 r is encoded 
by a sequence of pairs of finite trees (Lj, Ri)i^z- Clearly, this is an exchange- 
able sequence. Let atj be a permutation of the i-th and j'-th pairs; then 
cTjjT is a tree-parametrization of some other infinite Lorentzian triangula- 
tion (JijT, which only differs from T by a finite subgraph. By Lemma 15.11 
/3c(T) = (3c{o'ijT), and Hewitt-Savage zero-one law applies to /3c, considered 
as a function of the sequence {Li,Ri)i^z- Thus PdT) is constant a.s. 
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